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C^' Abstract 

^ I We discuss the qualitative features of the recent data on multiparticle production ob- 

served at the LHC. The tolerable agreement with Monte Carlos based on LO DGLAP 
^ ' evolution indicates that there is no qualitative difference between 'hard' and 'soft' inter- 

^ ' actions; and that a perturbative QCD approach may be extended into the soft domain. 

■ However, in order to describe the data, these Monte Carlos need an additional infrared 

5^ ■ cutoff /cmin with a value femin ~ 2 — 3 GeV which is not small, and which increases with 

collider energy. Here we explain the physical origin of the large /cmin- Using an alternative 
f"*) , model which matches the 'soft' high-energy hadron interactions smoothly on to pertur- 

bative QCD at small x, we demonstrate that this effective cutoff fcmin is actually due to 
the strong absorption of low kt partons. The model embodies the main features of the 
BFKL approach, including the diffusion in transverse momenta, In/cj, and an intercept 
^ ' consistent with resummed next-to-leading log corrections. Moreover, the model uses a 

. two-channel eikonal framework, and includes the contributions from the multi-Pomeron 

exchange diagrams, both non-enhanced and enhanced. The values of a small number of 
physically-motivated parameters are chosen to reproduce the available total, elastic and 
proton dissociation cross section (pre-LHC) data. Predictions are made for the LHC, and 
the relevance to ultra- high-energy cosmic rays is briefly discussed. The low x inclusive 
integrated gluon PDF, and the diffractive gluon PDF, are calculated in this framework, 
using the parameters which describe the high-energy pp and pp ^soff data. Comparison 
with the PDFs obtained from the global parton analyses of deep inelastic and related hard 
scattering data, and from diffractive deep inelastic data looks encouraging. 



1 Introduction 



The general-purpose Monte Carlo generators [1], like PYTHIA or HERWIG, describe the in- 
clusive spectra observed in hadron collider experiments assuming that there is some relatively 
hard parton-parton interaction in the central region, supplemented by the secondaries produced 
in backward evolution from the hard matrix clement to the incoming protons. It turns out that 
in order to reproduce the data, the infrared cutoff, fcmin, in hard matrix element should increase 
with collider energy reaching a value of about 3 GeV for = 7 TeV (as compared to 2.15 GeV 
at the Tevatron energy)^. Such a large value of k^in should be explained within a perturbative 
QCD framework. 

The most natural possibility is to say that in the low x region, relevant for high-energy 
collisions, the probability of rescattering, and/or additional secondary interactions, becomes 
large on account of the high parton density. Then the corresponding absorptive corrections 
(which are driven by a cross section a^^^ oc 1/kf) suppress the low contribution. 

Here we study this possibility at a quantitative level. We consider a model based on the 
conventional perturbative QCD in the high kt domain, which includes the main features of the 
BFKL approach [3] in a simplified form. In this way we describe the structure of the 'bare' 
Pomeron. Note that this QCD Pomeron has an internal variable: the transverse momentum, 
kt, of the intermediate partons. In order to extend the description into the lower kt domain 
we include multi-Pomeron contributions written in terms of Rcggcon Field Theory (RFT) [4]. 
Recall that RFT offers the description of high-energy 'soft' hadron-hadron interactions [5]. 
Thus we achieve a smooth transition from the perturbative QCD to the 'soft' domain. That is, 
we construct a model which describes both pure 'soft' and 'hard' phenomena within a common 
unified framework and with the same set of parameters. In other words, the model describes 
all facets of high-energy hadron interactions on the same footing. 



2 LHC data versus Monte Carlo predictions 

The inclusive single particle {ptiV) distributions measured at the LHC are in broad agreement 
with the Monte Carlo (MC) predictions, but there are significant discrepancies. In the most 
popular MCs, hke PYTHIA, even the soft interaction is simulated by a 'hard' subprocess (say, 
99 99:--) with 

kt > kmin = ko{s/soy, (1) 

accompanied by DGLAP evolution up to the scale, kt, corresponding to the hard subprocess; 
plus hadronization at the final stage. The parametric form of k^i^ is tuned to the data. Leading 
order (LO) evolution and LO PDFs are used. However, since the model accounts for exact 
energy-momentum conservation, the results include some NLO effects. 

^The numbers correspond to PYTHIA 8.1 [2]. 
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Figure 1: The theoretical formalisms appropriate to the various domains are indicated. As long 
as X is not too small, we have a well justified theory based on the perturbative QCD and DGLAP 
evolution. For smaller x we need to resum the terms enhanced by the large value of Inl/x inside 
DGLAP evolution, and to account for the semi-enhanced absorptive corrections generated by the 
Balitsky-Kovchegov (BK) equation [6]. Above the 'BK' curve more complicated multi-Pomeron 
diagrams enter and we cannot justify the result obtained by summation of 'fan' diagrams only, that 
is, the results based on the BK-equation. However, here we extend the partonic ladder structure of 
the Pomeron (generated by BFKL-like evolution in rapidity) to allow for a full set of multi-Pomeron 
exchange diagrams, and obtain a model which is applicable in the low kt region bounded by the 
dashed curve. Moreover, in principle, it is possible to use a more precise evolution kernel, which 
accounts for both BFKL and DGLAP logarithms, to cover the whole kt region in our approach. 



We now turn to the discrepancies. In the low pt < 0.5 GeV region the high-energy data 
exceed the MC predictions by, for example, up to 40-50% at pt — 100 MeV. To summarize, 
the data [7, 8] for the gi-integrated particle density dN/drj turn out to be larger, while the 
mean pt {{pt)) is smaller, than the predictions of the MCs. These observations indicate the 
presence of the configurations which violate the strong kt ordering, which is characteristic of LO 
DGLAP. Rather, at high energies, where we enter the low x domain, BFKL-like contributions 
are important, see Fig.l. That is, contributions which are ordered in x but not in kt- 

2.1 The 'old' Regge description of 'soft' data 

Before we introduce the BFKL amplitude, which sums up diagrams with so-called multi-Regge 
kinematics, we first recall the main features of the 'old' Regge description of soft high-energy 
hadron interactions, see, for example, [5, 9]. In this approach, the high-energy interaction 
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is described by Pomeron exchange, By 'Pomeron' we mean the rightmost singularity in the 
complex angular momentum or j plane. In the simplest case, the Pomeron is a pole in j-plane, 
which originates from the sum of ladder-type diagrams in which the transverse momenta of the 
particles are assumed to be limited, and not to grow with energy. 

Indeed, up to Tevatron energies the elastic amplitudes, are well described with a Pomeron 
which is a single pole at [10] 

af{t) = j = 1 + e + a't, with e = 0.08 and a' = 0.25 GeV~^ (2) 

The reason for the superscript 'effective' will become clear in a moment. On the other hand, in 
this one-pole approximation the single-particle inclusive cross section grows in the same way as 
the inelastic cross section, da/dydPpt — V5n(Pt)crinei, see [11]. In this case the particle density, 

dN _ 1 da 
dydpl (Jinei dydpl ' 
does not depend on energy, contrary to the data. 

The puzzle is solved by including the contribution of the multi- Pomeron exchange diagrams. 
Due to the cancellation^ it turns out that the single-particle inclusive cross section is still 
described by one Pomeron exchange, 

da/dydpl oc (4) 

whereas the growth of inelastic cross section is reduced by the contribution of the multi-Pomeron 
diagrams. In other words, e in (2), is an effective power^, which is smaller than A of the original 
'bare' Pomeron pole: 

a^^'\t) = 1 + ^ + a'p^'''\t). (5) 

In the present Monte Carlos these phenomena are included as the possibility of Multiple 
Interactions (MI) when a few pairs of incoming partons interact simultaneously, producing 
independent^ chains of parton showers. For each chain the evolution is convoluted with the 
corresponding 'hard' subprocess. 



2.2 Implications of the LHC data 

In the energy interval covered by the current LHC data, 0.9 < ^/s < 7 TeV, the particle density 
in the central region, = 0, grows as [7, 8] 

dN/dr] ~ (6) 

^This cancellation can be easily seen by applying the AGK cutting rules [12] to the (eikonal) multi-Pomeron 
diagrams. Eikonal diagrams correspond to Pomerons exchanged between the colliding protons. Later we will 
discuss so-called enhanced multi-Pomeron diagrams. 

^The 'effective' pole form, (2), only provides a reasonable description up to Tevatron energies. The elastic 
amplitude, Tei(s,6), already exceeds the black disc limit at small impact parameters for ^/s > 5 TeV, see, for 
example, [13]. 

■^Modulo the constraint of the conservation of the overall energy and momentum. 
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Taken together with the behaviour of the inelastic cross section (Jinei ~ s° °^, it means that the 
two-particle irreducible amplitude described by the 'bare' Pomeron pole, which according to 
AGK rules also specifies the inclusive cross section, leads to the result 

— = (Tinel^ S with A ~ 0.2. (7) 

dy dy 

Simultaneously, the mean transverse momentum of secondaries observed at the LHC [8] is 
found to have the behaviour 

{Pt) oc s'-'\ (8) 

In other words, the particle multiplicity, dN/dy, in the central region grows mainly due to the 
population of a larger phase space (oc (pt)'^) in the transverse-momentum distribution, while at 
very low pt the particle density dN/ dydpl is close to saturation. It should be emphasized that 
the energy dependence of mean (pt) is beyond description by 'old' Regge theory. 



2.3 Diffusion in In/c^ 

In Monte Carlos, based on LO DGLAP evolution, the transverse momenta of intermediate 
partons monotonically increase starting from the proton up to the matrix clement of 'hard' 
subprocess taking place somewhere in the central region. For the events with a really hard 
subprocess (characterised by a large scale) such configurations, ordered in kt-, indeed give the 
dominant contribution. However, in a normal inelastic event (with no large scale) there is no 
reason for ordering in kt. 

In particular, within BFKL evolution^ the transverse momentum can increase or decrease 
with equal probability, leading to diffusion in In/ct [15]. Thus, for the normal inelastic event, 
it would be natural to expect a large contribution from diagrams corresponding to 'partonic 
chains' not ordered in pt. 

It is surprising that the recent high-energy LHC data on inclusive single-particle distribu- 
tions are reasonably well described by the DGLAP-based Monte Carlos^. Partly, this may be 
explained by the number of parameters (which are not fixed by theory) used to tune the Monte 
Carlos. 

On the other hand, the fact that a rather large and increasing with energy infrared cutoff, 
(1), is required to describe the data, may indicate the suppression of the low kt contributions 
which enforces the increase of transverse momenta along the parton chain in going from the 
protons to the central plateau. 

^The evolution is started from some relatively large scale in order to justify pQCD and to neglect confinement. 
An example is an onium-onium interaction [14]. 

^Another surprise is that the reasonable description was obtained using the LO gluon PDF which grows with 
1/x much steeper than that in NLO case. Recall that experimentally the low x behaviour is fixed mainly by 
deep inelastic data where the heavy photon interacts with a quark. The absence of the 1/ z singularity in the LO 
quark-quark splitting function, Pqq{z), is compensated by a steeper 1/x behaviour of the LO gluon distribution. 
At first sight, the NLO PDF should be more relevant for the inclusive pp ^ X cross sections measured in the 
central region where the gluon-gluon interactions dominate. 
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Figure 2: Eikonal and enhanced multi-Pomeron contributions 
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Figure 3: The triple-Pomeron diagram, together with its ladder structure 



2.4 Low pt suppression and enhanced multi-Pomeron diagrams 

As mentioned in the Introduction, the suppression of low kt partons was actually expected. The 
absorptive cross section has the form a^^^ oc 1/kf. Therefore the absorptive effects stem the 
growth of 'wee' parton densities at low kf, while at larger kt the growth with 1/x is continued. 
In terms of Regge Field Theory (RFT) the absorption of intermediate partons is described by 
so-called 'enhanced' multi-Pomeron diagrams. 

So far, we have discussed only eikonal diagrams, with multi-Pomerons exchanged between 
the colliding protons. However, there is also the possibility that one or more Pomerons may 
couple to an intermediate parton. We then have a so-called (semi) enhanced multi-Pomeron 
diagram. The eikonal and enhanced diagrams are sketched in Fig. 2. The simplest enhanced 
diagram corresponds to the triple-Pomeron contribution shown in Fig. 3. Recall, that this 
diagram gives the leading contribution to the cross section for a proton dissociating into a 
high-mass system. The triple-Pomeron coupling^, g^p, is found [16, 17] to be rather small, 
93p/gN — 0.2, where is the Pomeron-proton coupling. Nevertheless, the final effect of the 
triple-Pomeron diagrams is rather strong at high energies, since the contribution of each diagram 
is enhanced by the large rapidity interval that is available for the position of the vertex. Indeed 
the enhanced diagrams lead to increasing amounts of high-mass dissociation with increasing 
collider energy. 



^(/3P = gf — g\ in the notation introduced below. 
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The simplest models contain only multi-Pomeron enhanced diagrams generated by triple- 
Pomeron vertices. On the other hand, such models lead to cross sections which asymptotically 
decrease with energy [18] . We note that there is no reason to assume that all the more compli- 
cated multi-Pomeron vertices, gl^ with n + m > 3, are exactly zero. In fact it is more reasonable 
to include these contributions. However, the values of the vertices, (7^, which couple n to m 
Pomerons, are not known forn-|-m > 3. In this paper we study the two most physically reason- 
able hypotheses. We formulate the hypotheses about the g^!^ behaviour in terms of conventional 
soft RFT in which the transverse momenta of partons are limited and where, for the moment, 
we omit the variable kt- 

The first hypothesis [19] is that 

(a) g:^^A{nm/2){g^Xr+--'. (9) 



This corresponds to a suppression 



exp(-Afi) (10) 



of each intermediate parton; like the absorption of the incoming beam particle in a nuclear 
target with the opacity Q. The second hypothesis is that 

(b) g:^^Gsp{gj,Xr+^-' with G3P = A, (11) 

where the last condition ensures that we are in the 'critical' Pomeron regime^. After the 
enhanced screening the 'renormalized' intercept (i.e. the position of the singularity in the j- 
plane) of the Pomeron is q;''''"(0) = 1. In this scenario the suppression factor at each intermediate 
parton vertex takes the form 

1 - exp(-AO) 

instead of the (10). 

For a smaller value of G^p < A the screened intercept a^^^{0) > 1. That is enhanced 
absorption is not sufficient to stop the power growth of the cross section. Only the eikonal 
diagrams restore unitarity. Actually this, 'supercritical' regime is like an eikonal model with 
the internal structure of the Pomeron slightly modified. On the other hand, a larger value of 
G3P leads to a decreasing cross section which dies out with energy. We regard these last two 
regimes as unrealistic, and do not consider them further. That is we focus on the 'critical' 
regime which is provided automatically in case (a), or with a specific value of Gsp in case (b). 
Further discussion is given in Section 5. 

Next, we would like to trace the smooth transition from the 'soft' low pt domain to the 
'semi-hard' region of relatively large pt- 



^More details are given in Section 5. 
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3 Matching the 'soft' and 'hard' contributions 



There is evidence of a smooth transition from the 'soft' Pomeron which describes 'soft' data 
to the perturbative QCD Pomeron describing the semi-hard region of relatively large pt, see 
subsection 3.1. This opens the way to link the description of high-energy soft interactions to 
the perturbative very low x, pt ~ few GeV domain, a region heavily populated by LHC data, 
see Fig. 1. We sketch how this is done in subsection 3.2 and give more details in the Appendix. 

3.1 Smooth transition between the hard and soft regimes 

There are phenomenological hints that at large distances the "soft" Pomeron should have qual- 
itatively similar structure as the "hard" (QCD) Pomeron, see [20]. Indeed, first, no irregularity 
is observed in the HERA data in the transition region, 

Q2 ~ 0.3 - 2 GeV^ between the 'soft' 
and 'hard' interaction domains; the data are smooth throughout this region. In particular, 
the observed cross sections for vector meson production, 7p — > Vp, are consistent with a 'soft' 
Pomeron of intercept ap(0) ~ 1.1 at low {Q^ + ^v)-, rising smoothly to Q;p(0) ~ 1.3 at large 
{Q^ -\- My). Next, a small slope a'p < 0.05 GeV~^ of the bare Pomeron trajectory, is obtained 
in global analyses of all available soft high-energy data, after accounting for absorptive correc- 
tions and secondary Reggeon contributions [21, 22]. This indicates that the typical values of 
kt inside the Pomeron amplitude are relatively large (a' oc 1/Af). Finally, recent 'soft' model 
data analyses [21, 22] which account for the enhanced absorptive effects find an intercept of 
the initial, bare Pomeron A = ap(0) — 1 ^ 0.3 close to the intercept of the BFKL Pomeron 
after the NLL corrections are resummed [23, 24]. Thus it looks reasonable to assume that in 
the soft domain we deal with the same perturbative QCD Pomeron; at least, there is a smooth 
transition from the soft to the hard Pomeron. 

The net effect is that the hare hard Pomeron, with a trajectory with intercept A = ap(0) — 
1 ~ 0.3 and small slope a', is subject to increasing absorptive effects as we go to smaller kt 
which allow it to smoothly match on to the attributes of the soft Pomeron. In the limited 
energy interval up to the Tevatron energy, some of these attributes (specifically those related 
to the elastic amplitude) can be mimicked or approximated by an effective Pomeron pole with 
trajectory intercept Aeff = (0) - 1 ~ 0.08 and slope a'^^ = 0-25 GeV"^ [10]. 

3.2 Outline of the procedure to hnk the hard and soft regimes 

We start with the partonic ladder structure of the Pomeron, F{ii, kt, b), generated by BFKL-like 
evolution in rapidity, see Fig. 4 



At each step of the evolution kt and the impact parameter, b, can be changed. It is important 
to note that, in comparison with RFT, we now have an extra variable, kt, that is the transverse 



dF{y,kt,b) 
dy 




(13) 
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Figure 4: The ladder structure of the bare Pomeron, F{y,kt,b), together with a symbolic sketch 
of the hadronisation process, cf (55). This figure shows the symbols for the various transverse 
momenta used in the text. 



momentum of the intermediate parton, in addition to the usual two 'soft' Regge variables: 
y = ln(l/x) and the impact parameter, b, which is conjugated to the transverse momentum, 
Qt, transferred through the entire ladder. 

Usually the evolution equation in perturbative QCD is written in terms of the amplitude 
integrated over the impact parameter b. However here we have included the explicit b depen- 
dence, since then it is easier to calculate the effect of the multi-Pomeron contributions. In 
general the kernel K depends on the difference b — b'. However, the b dependence is propor- 
tional to the slope a' of the bare Pomeron trajectory. Asymptotically the BFKL approach 
predicts a' — )■ 0. Indeed, analyses of soft data find it to be very small, a' < 0.05 GeV^, see 
[21, 22]. 

. For simplicity, we neglect this variation. Then the only b dependence of F comes from 
the starting distribution of the evolution, and not from the b dependence of K. 

The sum of the non-enhanced eikonal diagrams may be written in terms of the opacity fl. 
At each value of the impact parameter the result for the amplitude is 

T(6) = 1 -exp(-fi(b)/2). (14) 

Since at high energies the proton opacity Q is proportional to the gluon density, F{y,kt,b), 
unintegrated over the impact parameter b it is convenient to write the linear evolution equation, 
(13), in terms of F. Note that F depends not only on b and y (as in soft RFT), or only on y 
and kt (as in pQCD), but on all three variables. 

To include rescattering of the intermediate partons with the beam i and target k diffractive 
eigenstates, as in Fig 5, we have to include in (13) the absorptive factors (10) or (12), depending 
on one or another version of the model and to solve two coupled evolution equations. One 
evolving up from the target k aX y = 0, and one evolving down^ from the beam i ai y' = 

^The value of Yfc accounts for the fact that at larger k'f a smaller rapidity interval is available for the evolution. 
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n-y = 0withn = ln(s/A;;2), 



dy J Tik^ 

dF,{y') [£k[ 



dy' J nkf 



'2 



exp(-A[Q,(y) + Q,(y')]/2) K{kt,K) Fkiv)- (15) 



exp(-A[a(|/') + ^k{y)]/2) K{kt, k[) F,{y'), (16) 



where, for clarity, we have suppressed the fcj labels of the F's and Q's. Here we show the more 
familiar absorptive factors, exp(— Af2/2), of the amplitude at the intermediate parton vertices^*^ 
corresponding to hypothesis (a) of (9): in the case (b) it should be replaced by the form (12). 
The coupled evolution equations may be solved iteratively to give F{y,kt,b), for a whole range 
of fixed values of b. 

The opacities which enter (15,16) are given by the LO expression^^ 

Xnk{y,kt,b)= / Ny-;^asik,')F,{y,K,b), (17) 

where the integral over k^ accounts for the possibility of screening an intermediate parton c 
with transverse momentum kt by any additional ladder of smaller size^^, that is with k'^ > kf, 
see [25]. 

To allow for low- mass dissociation, p — >■ A?"*, we follow Good- Walker [26] and introduce 
diffractive eigenstates, i,k, which are those linear combinations of p,N*, ... which diagonalise 
the T matrix and only undergo elastic-like scattering. Thus we have to consider the interaction 
of the beam state i with the target state k; correspondingly the opacities Q in (15,16) are 
marked by the indices i and k. Two channels suffice [19]. 

Recall that the kernel K provides the possibihty of evolution in both directions - the mo- 
mentum \kt\ may be larger or lower than with equal probability. On the other hand the 
absorptive factor exp(— Af2(y, kt)) suppresses the production of a low kt partons. This acts like 
an infrared cutoff, /cmin, and makes the evolution asymmetric, leading to a larger probability of 
cvohition with kt increasing; that is, in the usual DGLAP direction. 

Besides this, we need to introduce a boundary condition at low kt, say kt = qo, close to the 
confinement region. We assume that the parton density is zero for kt < qo- Moreover, to avoid 
possible double counting, and to fix the boundary between the low- and high-mass dissociation, 
we introduce a threshold Ay = 1.5 in rapidity at which we start the evolution of (15,16). That 



^"Since we are dealing with the amplitude F, and not with the cross section, we use here exp(— Af2/2) and 
not exp(— Ail). 

^^The normalisation of this form has been adjusted such that F{y,kt,b) is the doubly-unintegrated gluon 
distribution, as given in (25) below. Wc use the LO pQCD couplings with one-loop as (n/ = 3, A = 150 MeV) 
and allow for some renormalization due to NLO effects. 

^^At a small k'^ < kt the gluon with a large Compton wave interact coherently with the colour charge of the 
whole (colour neutral) ladder. This contribution is negligible. 
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Figure 5: A typical enhanced multi-Pomeron exchange diagram 



is, we start the upwards evolution at y = yo = Ay = 1.5, and not at ?/ = 0, and the downwards 
evolution from y' = yo = Ay = 1.5, and not from y' = 0. The relation between y and y' is 

y' = Yk~y where Yk = ln{s/k^), (18) 

recall kt > qq. y^ specifies the low- and high-mass dissociation regions, and high-mass dis- 
sociation takes place in a rapidity interval ln{s/kf) — 2 Ay. Low-mass dissociation is allowed 
for via the diffractive eigenstates, and high-mass dissociation is included via the enhanced 
multi-Pomeron exchange diagrams. 

3.3 The parameters of the model 

The idea is to study, in a semi- quantitative way, the main features of soft and semi-hard data in 
terms of a realistic model with just a few physically-motivated parameters, and not to perform 
a multi-parameter x^-analysis of the data. In this way we hope that we can provide a better 
understanding of the physics which underlies the description of the data. Using the procedure 
outlined above, we start from a formulation in the larger kt (perturbative QCD) domain, and 
then extend the description into the lower kf region in order to reproduce the 'soft' data. 

It is informative to gather together the small number of parameters, that we have introduced 
above (or have implied), that are to be tuned to describe all the gross features of the available 
'soft' high-energy data. There is basically one parameter (or sometimes two) that is mainly 
responsible for each phenomena: 

• A = ap{0) — 1 giving the intercept of the bare Pomeron trajectory, which controls the 
energy dependence. In other words, the kernel K is normalized to give a final intercept 
A for the bare Pomeron before absorptive effects are included. The slope a'p of the 
trajectory is set to zero, since it is found to be very small. 

• d which provides a smooth transition between the logarithmic contributions with kt <^ k't 
and kt ^ k't. It thus controls the BFKL-like diffusion in InA;^. 

• A^' which specifies the initial gluon density, and in this way fixes the Pomeron-proton 
coupling. 
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• the parameters Ci and C2 which determine the proton radius and the shape of the corre- 
sponding form factor. 

• A ( or rather f^p of (19) below) which determines the strength of the triple- (and multi-) 
Pomeron couplings, which are constrained by data on high-mass diffractive dissociation. 

• 7j which specify the Good- Walker diffractive eigenstates [26], and which are determined 
by low-mass diffractive dissociation (the coupling of state i to the Pomeron is proportional 
to 7i). 

• Ho which, in order to avoid double counting, separates low- and high-mass diffraction (the 
low mass dissociation is written in terms of the Good- Walker eigenstates while the high 
mass dissociaton is described in terms of multi- Pomeron (RFT) contributions). 

• qo, the infrared cutoff, which together with N', controls the absolute value of the bare 
one-Pomeron exchange cross section. 

Except for the important introduction of an explicit treatment of the kt dependence, all 
other features of the model are practically the same as described in our previous publications 
[19, 21]. 

3.4 Simultaneous description of soft and semi-hard data 

Having tuned the parameters to describe the available high-energy 'soft' data {atot, da^i/dt, 
dasB / dtdM^ , ..), we check whether the model can simultaneously describe semi-hard phenom- 
ena, such as the inclusive single-particle pt distribution, the gluon PDF and diffractive gluon 
PDF in the low x domain. Before presenting the detailed results, it is informative to first gain 
some insight into how the model offers the opportunity to simultaneously describe such diverse 
data. 

3.4.1 Choice between hypotheses (a) and (b) for multi-Pomeron coupUngs, g'^ 

We saw, in subsection 3.2, that the amplitude, F{y, kt, b), is obtained from a perturbative QCD 
Leading Logarithmic (LL) formalism with an infrared cutoff, qq. The low kt region, kt ^ qo, 
gives the 'soft' contribution, while at larger kt we have the conventional perturbative QCD 
result. 

Note that in our approach the same parameter A simultaneously describes (i) the absorptive 
effects which suppress the low kt contributions, so pushing the majority of partons into a larger 
kt region, and (ii) proton dissociation into high-mass systems with ^ m\ since the value 
of A determines the size of the triple- and multi-Pomeron contributions. Data corresponding 

to these processes, da^D / dtdM"^ , are available at Tevatron [27, 28] and lower energies [29]. As 
noted above, these are the data which fix the value of A, or rather f^p defined below. In this 
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way we obtain a parameter-free prediction for low kt suppression and thus for the whole kt 
behaviour (modulo simplifications used in the present model). 

Clearly, the dependence (the 'additional' variable) along the ladder plays an important 
physical role. Note that now A cannot be considered just as a constant number. Its value 
depends on kt- It is reasonable to use the LL expression corresponding to the BFKL triple- 
Pomeron vertex 

A = fspNcasihy , (19) 

where the constant /sp, which allows for the possible renormalization caused by the next-to- 
leading corrections, may be regarded as the new parameter in place of A. Clearly we expect 
fsp ~ 0(1). A is dimensionless, since the absorption during the evolution is written in terms 
of the opacity fl.{y, kt, b). The small probability of large kt partons rescattering is due to the 
decrease of Q{kt) with increasing kt. 

Analogously, the starting distribution for the evolution of F is written in terms of the LO 
(Born) cross section, which decreases as l/k^, multiplied by the 'input' parton density 

Fiyo,klb) = N'mNcOCs{kt)/T:kl . (20) 

Here N' is the normalization corresponding to the number of input gluons at a; = exp(— t/o), 
and (5{h) describes the parton distribution in the impact parameter plane. It is given by Fourier 
transform of the effective two-gluon form factor /3(t) 

m = ^ / e^^^-^md'Qt. (21) 

where t — —Qt, and where we use the parametric form 

/3(t) = e-V(l-t/ci)2. (22) 

Recall that, due to colour coherence, a parton with high kt in the ladder cannot be screened 
by one with lower kf. The lower kf amplitude, with large wavelength, 'sees' a point-like colour- 
less object at high kt, and does not interact with it. Therefore the 'effective' value of the 
triple-Pomeron vertex, and the strength of suppression (due to the multi-Pomeron interactions) 
decreases with increasing kt- 

The major contribution to the total and diffractive dissociation cross sections comes from 
the lowest kt, close to the infrared cutoff qq. In the case of crtot it is possible to compensate a 
larger value of go by a larger input parton density N'. However the relative value of high- mass 
single dissociation, dasD/dM"^ still decreases if the value of go is increased. 

Now, the multi-Pomeron couplings specified by hypothesis (a), given in sTibsection 2.4, 
correspond to strong enhanced absorptive factors of the form cxp(— Af2). For this reason, we 
need to choose a rather low infrared cutoff go = 0.5 GcV to obtain a satisfactory description 
of (Tel, daei/dt, and high-mass dissociation, dasn/dM^. Since the normalization is fixed by 
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the soft cross sections at low pt ~ a-nd since all the higher contributions are suppressed 
by the ratio ql/Pf, we fail to reproduce the inclusive single-particle spectra, da/dydp^, in the 
larger pt domain^^. Indeed, the model underestimates the inclusive cross sections measured at 
CERN and the Tevatron by about factor of 3-8 for pt > 4: GeV. Simultaneously, the diffractive 
gluon PDF predicted by the model turns out to be a factor 3-5 smaller than that determined 
by the diffractive deep inelastic data measured at HERA. These are rather general features of 
hypothesis (a). They do not change too much under variation of the parameters of the model, 
including the structure of the Good- Walker eigenstates [26] used to reproduce low-mass proton 
dissociation. Therefore below we will concentrate on the case (b). 

With the multi-Pomeron couplings specified by hypothesis (b), the suppression of the low 
Pt domain is weaker. Here, a good description of 'soft' cross sections, (Jtot, daei/dt, dasn/dM^, 
is obtained with the infrared cutoff choscn^^ to be qq = 0.85 GcV. This increase of the cutoff 
from 0.5 to 0.85 GeV also makes the description of the single-particle inclusive cross section 
much better, see Fig 6. Now the predicted values da/dydpl are consistent with the present 
collider data for ^ 5 GeV. At lower pt the naive prediction, based on a simple fragmentation 
of the gluon minijets, without a more precise treatment of hadronization, still underestimates 
the single-particle inclusive cross section. However, the problem of the deficit of low pt hadrons 
may be solved by accounting for the particles produced via hadronization (that is the breaking) 
of the colour (LUND) strings. The role of this effect cannot be evaluated analytically. Here, 
we need a Monte Carlo. 



3.4.2 The gluon PDFs: g(x,iJ,'^) and g^(xp, z, /j,'^) 

Another check of how well our model can be matched to perturbative QCD, are the predictions 
it gives for the gluon PDFs. The 'effective' gluon distribution may be calculated assuming that 
each kt component corresponds to a ladder with a gluon of transverse momentum fcj. Such a 
ladder gives us the gluon distribution F{y, kt, b) unintegrated over kt and b. Indeed, in terms 
of the unintegrated gluon distribution fg(x,k^, ii"^), the integrated distribution is 

dk^ 

xg{x,ij,^) = J -^fg{x,klii^) , (23) 

where the gluon density, fg, unintegrated over kt collects the gluons from the whole impact 
parameter b plane. That is 



fg{X, kllJ?) 

and 



= jd'bF{x,klb,i,') (24) 



xg{x,i^^) = Jd'bJ' dklF{y,kt,b), (25) 



^^The formula for the calculation of the pt distribution are given in (54) of Appendix B.3 
^^This value of the infrared cutoff was used in a series of Durham papers for exclusive Higgs boson production 
and other Central Exclusive Diffractive processes [31]. 
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P, (GeV) 

Figure 6: The description of the Tevatron data [30] for the single-particle pt distribution, and the 
prediction for the LHC energy of 14 TeV. The continuous and dashed curves correspond to choices 
(ii) and (i), respectively of the diffractive eigenstates, see Section 4.1. 



where y = Inl/x. However, care needs to be taken with this identification. In order to compare 
with the gluons obtained from global parton analyses, that is gluons obtained from fits to DIS 
and related 'hard' data^^, we have to calculate the unintegrated distribution fg in the target 
proton neglecting the absorption caused by the beam opacity fli. Thus we must use an equation 
of the form of (15) but with fli = 0, such that there is no absorption caused by the beam proton. 
We denote the target opacity obtained in this way by fl, and the unintegrated gluon density 
by F. Hence the integrated gluon density is actually 

xg{x,ii^) = J]|afc|2 f (fb r dklFiy,kt,b), (26) 

k 

where the sum over the Good- Walker eigenstates k in the incoming proton wave function 
\p) — Tiak\k) is included^^. 

The so-called diffractive gluon distribution, g^{xp, z, jj?'), can also be predicted. In order to 
calculate g^{xp, z, //^) we have to replace the full opacity in (26) by the opacity corresponding 
to events with a rapidity gap between the target k and an intermediate parton c placed at 
y = yp = ln(l/a;p), that is, to events with elastic c — k scattering. For the favoured hypothesis, 
(11), of the behaviour of the multi-Pomeron couplings, gl^, the starting value of at y — yp 

^^We emphasize that the parameters of our model were tuned to 'soft' pp and pp data, not to DIS data. 
^^AU the coefficients in the decomposition \p) = ^ ak\(t>k) are taken to be Ofe = l/v^ (with k = 1, 2). 
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is given by 

F,^(yp, k„ b) = y (1 - e-^^^(^-'=-^)/^)' , (27) 

while for hypothesis (9) 

F^{yp, kt, h) = Ae-^^*(^^''=*'^)/' (l - e-^^*(^^''=*'^)/2) ^^(yp, fc*, 6)/2 . (28) 

After this, the diffractive density F-^ is evolved according to (15) with VLi = 0. 

The predictions of the gluon PDFs using the model, with the favoured hypothesis (b) for 
the multi-Pomeron couplings, are also encouraging for the attempt to obtain a simultaneous 
description of soft and semi-hard phenomena. The diffractive gluons, g^{xp, z, fi^), predicted 
by the model are in broad agreement with those obtained from the NLO analyses of the HERA 
diffractive deep inelastic data. Since the model accounts for enhanced absorptive corrections 
(which is a non-linear effect), the diffractive gluons grow with the scale //^ a bit more slowly 
than those obtained by a NLO analysis based on linear DGLAP evolution. The same is true 
for the inclusive integrated gluons g{x,ii^). Actually the inclusive gluon distribution predicted 
by the model turns out to be between the values obtained from the LO and NLO global parton 
analyses of deep inelastic and related hard scattering data. As discussed above, it is surprising 
that the spectra of secondaries are better described by the Monte Carlo models which use LO 
and not NLO gluons. The detailed comparison of the predictions of the gluon PDFs with the 
previously known values is given in Tables 3 and 4 of the next Section. 

To conclude this Section, we emphasize that both the global inclusive and the diffrac- 
tive gluon distributions were not fitted, but are calculated in the model, based on the val- 
ues of the parameters which provided the best description of the purely soft cross sections 
((7tot, daei/dt, dasu/dM'^). That is, within the model, the same values of the parameters con- 
trol the X- and scale-behaviour of the gluon PDFs, as well as the suppression of the inclusive 
single-particle distribution at lowpt, as discussed in subsection 3.4.1. 

4 Results of the fit to soft data and predictions for semi- 
hard phenomena 

In this section we present the description of the 'soft' data obtained by tuning the parameters 
of the model, as well as the detailed predictions for the gluon PDFs. The soft data on their own 
do not definitively distinguish between hypotheses (a) and (b) of Section 2.4 for the form of the 
multi-Pomeron couplings. However, as explained in the previous Subsection 3.4, the two sets 
of predictions for the semi-hard phenomena strongly favour hypothesis (b), with the infrared 
cutoff qo — 0.85 GeV. We therefore present results for just this hypothesis. First, we give the 
values of the parameters that we obtain. 
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4.1 The parameters and their values 



The list of the, remarkably few, parameters of the model was given in Section 3.3. Moreover, 
the parameters are physically motivated, and for many we have a good idea what their value 
should be. We have already defined A in (5), N' in (20), Ci and C2 in (22), A (which was 
replaced by fsp) in (19), qo and yo at the end of subsection 3.2. The diffusion parameter d will 
be specified in (37) below. 

Now let us define the parameters 71 and 72, which are needed to account for low-mass proton 
dissociation. The incoming proton wave function is written as a sum of two components^'': the 
Good- Walker diffractive eigenstates, i = 1,2. These states are preserved (and not mixed) by 
an interaction with the Pomeron. That is, they undergo elastic-like scattering. The ratio of 
the Pomeron couplings to the first and the second state was chosen to reproduce the observed 
value of the low-mass single dissociation cross section^^ crs°D ^ — 2 mb at the CERN-ISR energy 
of ^/s — 31 GeV [32]. However, still some freedom remains. In general, each eigenstate may 
have its own form factor, which specifies the distribution of the partons in b space. We assume 
that the distribution has the same form for each state, except for the radius of the state, Ri. 
We consider two possibilities: 

(i) i?f oc (Tj, that is, the same parton density at the origin {b = 0) of each state, 

(ii) Ri oc (Tj, which corresponds to a BFKL-like coupling f3i oc l/fcj. 

In both cases, after tuning the parameters, we obtain qualitatively (and mainly quantita- 
tively) the same results. The main difference is that the model with the BFKL-type couplings 
of (ii) predicts the diffractive PDF of the gluon to be closer to the conventional values obtained 
from the diffractive HERA data. Although (ii) is the favoured choice we present the results for 
both possibihties for the 7i's. The values of the parameters are hsted in Table 1. 



4.2 The description of the 'soft' data 

The values of the total, elastic and single proton dissociation cross sections are given in Table 
2 for collider energies ^/s = 1.8, 7, 14, 100 TeV. The detailed formulae used to calculate the 
cross sections are given in Appendix B. 

We show in Fig. 7 the quality of the description of the data for the elastic differential cross 
section. In the tuning of the parameters to describe these data, we find that the slope of 
the ^-distribution at CERN-ISR energies is a bit too steep, and that it is a bit too flat at the 
Tevatron energy. This indicates the need to add a small non-zero slope, a', of the bare Pomeron 
trajectory, which was neglected in the present computations^^. We also show the prediction for 
differential elastic cross section at an LHC energy of ^ — 14 TeV. 

-'^^A two-channel eikonal was shown to suffice, see [19]. 

^^Hcrc, crg°Q ^ corresponds to data with the mass of the dissociating system M < 2.5 GeV, which approxi- 
mately equates to y < 1.5. These data are, to our knowledge, the only available measurement of CTg^^^, and 
even this has large uncertainties. Thus the identification of the diffractive eigenstates is far from unique. 

^^In particular, a small non-zero contribution to a' is generated by the pion-loop insertion into the Pomeron 
trajectory [33, 34]. 
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A 


u.o ± 


u.oo 


iV 


p. 8 
O.o 


O.oO 






Q 

o 




0.1 


0.9 




0.6 


1 


7i 


1.57 


1.5 


72 


0.43 


0.5 



Table 1: The values of the parameters for assumptions (i) and (ii) for the diffractive eigenstates. 
The other parameter values are — 2, yo — 1.5 and go = 0.85 GeV; also GeV units are used for ci 
and C2. 



energy 


Ctot 




^lowAf 
^SD 


highM 
^SD 


tot 
^SD 


_lowM 
^DD 


1.8 


72.8/72.5 


16.3/16.8 


4.4/5.2 


7.0/7.8 


11.4/13.0 


0.3/0.4 


7 


89.0/86.8 


21.9/21.6 


5.5/6.7 


9.9/10.2 


15.4/16.9 


0.5/0.7 


14 


98.3/94.6 


25.1/24.2 


6.1/7.5 


11.5/11.3 


17.6/18.8 


0.6/0.9 


100 


127.1/117.4 


35.2/31.8 


8.0/9.9 


16.7/14.4 


24.7/24.3 


0.9/1.6 



Table 2: Cross sections (in mb) versus collider energy (in TeV). The first number corresponds to 
choice (i), and the second number corresponds to choice (ii) of the diffractive eigenstates. usd 
denotes the sum of single-dissociation of the beam and the target. 

The description of the total cross section data is shown in Fig. 8(a). The screening correc- 
tions, which arise from the 'enhanced' multi-Pomeron diagrams, that is, from the high-mass 
dissociation, slow down the growth of the cross section with energy. Thus, the model predicts 
a relatively low total cross section at the LHC energy of 14 TeV: 

cxtot ^ 95 - 100 mb. (29) 

The energy behaviour of the elastic and diffractive cross sections are shown in Table 2 
and Fig. 8(c,d), where the high-mass single dissociation cross sections include all events with 
M^/s < 0.05. The dependence of the cross section for high-mass dissociation, ^(Pa/dtdC,, on 
^ = M'^/s is compared with the Tevatron CDF data [27, 28] in Fig. 8(b). Recall that in 
our model we have not included the secondary Reggeon contributions, and the corresponding 
results are shown by the dashed lines. The continuous lines include the secondary Reggeon 
contributions computed as in [21]. 

4.3 Proton single dissociation at the LHC 

The energy behaviour of diffractive cross sections is shown in Fig. 8(c,d) for the two choices 
of radii, (i,ii), of the Good- Walker diffractive eigenstates. Due to a small probability, S^, that 
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Figure 7: The t dependence of the elastic cross section, and the prediction for 14 TeV. The bolder 
and fainter (red) curves correspond to choices (i) and (ii) of the diffractive eigenstates respectively. 
The references for the data are the same as those given in [34]. 



the rapidity gap survives rescattering in the centre of the disk, the cross section of diffractive 
dissociation comes mainly from large impact parameters, b, that is from the periphery. Now, 
in the favoured model (ii), the component with a large inelastic cross section has a much larger 
radius and a smaller parton opacity. Therefore, model (ii) predicts noticeably larger diffractive 
dissociation. 

At this point, it is worth emphasizing again that the available soft data at pre-LHC energies 
are not sufficient to fix the details, and the precise values of the parameters, of the model. 
Another example of this, is the Ostapchenko model [35], which is rather close ideologically to 
our approach^°. It also fits the pre-LHC data. Nevertheless, with set (A), Ostapchenko predicts 
that (Jtot is about 30 mb larger, and cxgjf^*^ is more than twice smaller, than our expectations 
at ^/s = U TeV. 

As mentioned in Subsection 4.1, the identification of the diffractive (Good- Walker) eigen- 
^"The Ostapchenko model is discussed in Section 5. 
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Figure 8: The description of the data (a) for the total pp{pp) cross section, and (b) for single 
dissociation. The bolder and fainter (red) curves correspond to choices (i) and (ii) of the diffractive 
eigenstates respectively. The dashed lines in plot (b) correspond to the result if the secondary 
Reggeon contributions were to be neglected. Plots (c),(d) show the energy dependence of elastic 
and diffractive cross sections for choices (i),(ii). ctsd denotes the sum of single-dissociation of the 
beam and the target. 
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^ 5 

Figure 9: Cross sections for single proton dissociation, ^dasD/d^ integrated over t, for two LHC 
energies. The bolder and fainter (red) curves correspond to choices (i) and (ii) of the diffractive 
eigenstates respectively. The dashed curves are obtained if the contributions of the secondary 
Reggeons are omitted. The secondary Reggeon contributions are computed as in [21]. These 
contributions are not, of course, included in the values of ctsd listed in Table 2; formally, they are 
not of diffractive origin. 



states is a major uncertainty, since we have little experimental information on low-mass disso- 
ciation. Measurements of low-mass proton dissociation at the LHC can possibly be obtained 
using the Zero Degree Calorimeter and the Forward Shower Counters [36]. This would shed 
light on the structure of the diffractive eigenstates and help to better constrain the model. 

The expected cross sections of high-mass single diffractive dissociation, integrated over the 
transverse component of the momentum transferred to the proton (that is, integrated over t), 
are shown in Fig. 9 at the LHC energies of 7 and 14 TeV. Here we have plotted the results 
down to very small ^; ^ ~ 10^''. It is clear that measurements in the region ^ ^ 10~^ would 
be very valuable. At present such small ^ values cannot be measured by the forward proton 
detectors, but we may expect that such events with very large rapidity gaps can be selected, 
using forward calorimeters and, hopefully, the Forward Shower Counters [36]. For ^/s = 14 
TeV, ^ = 10~^ corresponds to the production of a diffractive system of the mass M = 14 GeV. 
The position of the edge of the large rapidity gap, y — 4, may be observed in the calorimeter. 

An interesting, but a little unusual, feature of the plots 9(a) and (b) are the maxima at 
^ ~ 3x 10~^ and ^ ~ 10~^, that is, at M ~ 12 and 14 GeV, for the two LHC energies respectively. 
In the simplified triple-Pomeron approach one expects a growth of M'^da/dM'^ oc 1/(M^)^ with 
M decreasing; we assume here a Pomeron intercept a{0) = 1 + A. Such behaviour (with the 
intercept 'renormalized' by more complicated, enhanced diagrams) is, indeed, observed in our 
model in the ^ = 10~^ — 10~^ interval. However, at very small ^, the probability of dissociation 
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//2 — 4 


//2 — 1 6 


,,2 _ aA 


1/2 — 256 

JL* ^ <J W 


T — in-2 

Ju — ±.\J 




9 8/10 1 


15 3/159 


99 1 /9'? 9 




(3.3-5) 


(5.2-6) 


(6 3-7) 


(7 1-7 5) 


X = 10-3 


6.7/6.9 


14.1/14.7 


24.4/26.0 


39.0/42.0 




(3.8-9) 


(9.6-16) 


(15-22) 


(19.6-27) 


X = 10-^ 


9.1/9.3 


21.0/22.1 


39.9/43.2 


68.8/75.8 




(4.2-16.5) 


(16.1-36) 


(29.8-56) 


(44.4-75) 



Table 3: The values of the integrated gluon distribution, xg{x,ii'^) with /i^ in units of GeV^, 
computed using the soft model, as given by eq.(26). As in Table 2, the two numbers correspond to 
choices (i) and (ii) for the diffractive eigenstates. For comparison the numbers in brackets correspond 
to the integrated gluon distribution determined in the MSTW NLO (first) and LO (second) global 
parton analysis [37]. 

starts to decrease with decreasing ^. The explanation is that the cross section coming from low 
kt components is strongly suppressed by enhanced absorptive effects. Thus most of the high- 
mass dissociation comes from relatively high kt. On the other hand, at small M, the rapidity 
(In 1/x) interval and available phase space is insufficiently large to generate high kt partons 
starting from the evolution of 'soft' partons in the proton wave function. This explanation 
is verified by using the same model, but without the kt dependence^^ . Then the maximum 
disappears and ^da/d^ continues to increase with decreasing ^. The growth of the cross section 
for high-mass dissociation with increasing M, at a relatively low M^, is a qualitatively new 
feature of a model which explicitly accounts for the kt dependence of the absorptive effects 
during the evolution. 

4.4 The gluon PDFs: g{x,fi^) and g^{xp, z, /i'^) 

In Tables 3 and 4 we present the predictions for inclusive and diffractive gluon densities respec- 
tively, obtained from the parameter values listed in Table 1. A discussion of the comparison 
of the predicted values of the gluon PDFs with the conventional determinations was given in 
subsection 3.4.2 

4.5 Rapidity gap survival probability 

The above model also allows us to calculate the probability that the rapidity gaps in diffractive 
processes survive both eikonal and enhanced rescattering. Recall that the absorption caused 
by rescattering of intermediate partons is included as the factors (12), or (10), in equations 
(15,16) for the opacities, depending on whether hypothesis (b), or (a), for gl^ is adopted. So we 

Allowing only a very small interval of kt variation. 
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fO 24) 


(0.53) 


fO 72) 


(0.85) 


Xp = 
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0.24/0.40 


0.33/0.56 


0.40/0.70 


0.49/0.86 


z = 


0.2 


(0.31) 


(0.54) 


(0.65) 


(0.68) 


Xp = 


0.001 


0.27/0.44 


0.40/0.69 


0.56/0.96 


0.75/1.32 


z — 


0.05 


(0.32) 


(0.76) 


(1.09) 


(1.30) 



Table 4: The values of the diffractive gluon distribution, xpzg{xp, z, /x^) with /x^ in units of GeV^, 
computed using the soft model, as given by eqs.(26,28). As in Table 2, the two numbers correspond 
to choices (i) and (ii) for the diffractive eigenstates. For comparison the numbers in brackets 
correspond to the diffractive gluon distribution determined from the MRW analysis of the diffractive 
HERA data [38]. 

have the possibihty to quantify the role of the suppression of processes with rapidity gaps due 
to enhanced rescattering, by solving the evolution equations with and without the enhanced 
absorption factors. We denote the 'target' k opacity, calculated without the absorptive factor 
from the beam i side, as 1^^; and the corresponding matrix element as M_. Further, if we denote 
the matrix element which includes the absorptive factors as M.^^^, then we obtain [39] 

'S'enh = I ^ |2 • (30) 

Note that now the value of S^^^ can be calculated for any fixed virtuality (transverse momenta) 
of the partons which initiate the hard subprocess. Then it can be included inside the kt integral 
for the hard subprocess amplitude. In this way, we can account precisely for the 'soft-hard 
factorization breaking' in central exclusive production and other processes. 



energy 


(i) 


(ii) 


7 TeV 
14 TeV 


0.013 
0.008 


0.024 
0.015 



Table 5: The effective (corresponding to the slope B — A GeV~^) survival factor, {S"^), of the 
rapidity gaps for the exclusive production, pp ^ p + H + p, of a Standard Model Higgs boson of 
mass 120 GeV at two LHC energies. Both eikonal and enhanced rescattering are taken into account. 

Indeed, our model allows the calculation of the rapidity gap survival factor, S^{b), for any 
diffractive process as a function of the impact parameter b. A topical, and important, example 
is the central exclusive production of a Higgs boson via the process pp p + H + p, where 
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the '+' signs denote large rapidity gaps. The overall survival factors (-S^), are given in Table 
5 at two LHC energies , for two choices, (i) and (ii), of the diffractive eigenstates; here the 
survival factors take into account both eikonal and enhanced rescattering. Recall that (ii) is 
the favoured choice. The values are in agreement with our previous determination at ^/s—14 
TeV [39] 

{S') = OM5t'o'oL- (31) 

Note that, in order to compare with other results^^, we have presented the effective values of 
{S"^) which correspond to the t-slope B = A GeV^^ of the bare hard cross section, that is, 
before accounting for absorptive effects. The value of {S"^) increases with the slope B, but the 
ratio {S"^) / B"^ is approximately constant for reasonable variations of B] again more discussion 
is given in [39]. 



5 Comparison with the other approaches 

At the moment, the only other approach which includes a complete set of the multi-Pomeron 
vertices is the model presented by Ostapchenko [35]. A comparison with our model is instruc- 
tive. As in our case, the Ostapchenko model uses a 2-channel (eikonal) Good- Walker formalism 
to account for low-mass diffractive dissociation, and the approach takes some account of the 
internal transverse-momentum structure of the Pomeron. 

First, recall, that in perturbative QCD, the BFKL Pomeron is not a pole, but a cut, in the 
complex angular momentum plane. This reflects the non-trivial structure of the BFKL vacuum 
singularity. Such a Pomeron may be treated as a series of components of different transverse 
size. Ostapchenko models this feature by simply introducing two different poles corresponding 
to small- and large-size Pomerons, whereas here we allow for the explicit kt dependence of the 
wave function corresponding to the t-channel Pomeron state; that is, we allow for the BFKL 
cut. In other words, the two poles of the Ostapchenko model may be considered as a way 
to mimic the BFKL cut, while we study the full kt dependence. Moreover in [35], the two 
poles were combined into a single propagator. As a result, one cannot consider the effect of 
a single Pomeron pole, with the contribution of the other component decoupled. This is a 
shortcoming of the model. In particular, in multi-Pomeron diagrams, there is screening of the 
small-size Pomeron by the large-size Pomeron, while in perturbative QCD this effect is absent. 
To be more precise, at leading order, the triple- Pomeron vertex contains (after the azimuthal 
integration) a 0-function (see [25]), which reflects the fact that a long- wavelength gluon (from 
the large-size Pomeron) interacts coherently with the small-size dipole formed from a pair of 
t-channel gluons in the ladder which represents the small-size BFKL Pomeron. Therefore the 
coupling is proportional to the whole colour charge of this small-size dipole. Since the Pomeron 
whole is colourless, the vertex vanishes. 

^^Notc that the very recent vakie of the survival factor obtained by the Tel- Aviv group [40] is an order-of- 
magnitude greater than their previous determination [22], and more in agreement with our value. 
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Contrary to our previous soft models [21, 34], we have seen above that the introduction 
of constraints from semi-hard phenomena now favours multi-Pomeron couphngs g^^ given by 
hypothesis (b) of (11), which, for the convenience of this discussion, we write as 

9l = r,p{g^\r^^-\ (32) 

with = (yfAfAA. The advantage of form (11,32) is the possibihty to use the simplest gener- 
alisation of the AGK cutting rules [12]. These rules determine the relative sizes of the contri- 
butions of the processes which result from the different cuts of the exchanged Pomerons (in a 
multi-Pomeron diagram) in terms of combinatorial factors only. 

This same form (32) was used in the Ostapchenko and KPT [41] models. However Ostapchenko 
uses two parameters: r^p for the triple-Pomeron vertex and A to allow for the other multi- 
Pomeron couplings, whereas we consider the so-called 'critical Pomeron' regime in which the 
value of Tzp — Qn^^ is fixed. 

An attractive feature of our critical Pomeron approach is that the two-particle irreducible 
amplitude, fl{b) of (14), generated by the evolution equations (15,16), never grows as a power 
of energy, but instead 

n{b)(x\ns. (33) 

On the contrary, if the general form (32) with an arbitrary value of the coupling r^p were to 
be used, then a power-like asymptotic behaviour, ft oc is obtained with a renormalized 

intercept (see [41] and Sect. 4 of the first paper of [35]) 

ar = «P " - rsp/^ivA (34) 

Here and ctp" denote the Pomeron intercepts before and after including the enhanced 

multi-Pomeron diagrams. If Op'^ < 1 then the cross section asymptotically decreases with 
energy. On the other hand, if the renormalized intercept «p" > 1 then unitarity is still not 
satisfied by the renormalized Pomeron exchange. In other words the enhanced contributions 
arising from interactions between the partons inside one parton cascade are not sufficient to 
provide saturation and to stop the growth of parton densities. Unitarity is only satisfied after 
eikonalization of the final amplitude. 

It is worth emphasizing that when parameters were tuned to describe the data, the Ostapchenko 
model resulted in an intercept ap^ < 1 for the large-size Pomeron component, while the inter- 
cept of the small-size component became very close to one, for example, for set C it was found 
that ap"^ — 1.03 [35]. This would indicate that Nature prefers the so-called 'critical Pomeron' 
regime (with ccp" = 1 resulting in a logarithmic energy behaviour of the irreducible amplitude), 
which is generated automatically by our form of the multi-Pomeron vertices. 

Note also that models which allow only triple-Pomeron vertices, such as that used by the 
Tel- Aviv group [22] , predict an asymptotically decreasing cross section. The reason is that (34) 
gives a negative ctp"^^ in the limit A ^ 0. However, the MPSI approximation used in [22] leads 
to a constant cross section, which indicates that the method is only applicable to a limited 
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energy interval, see also [35] . There should be the parameter which controls the validity of the 
approach of Ref. [22]. Besides this the Tel- Aviv model accounts for only one component of the 
Pomeron and its kf structure was not considered at all. 

6 Implications for ultra-high-energy cosmic rays 

A crucial ingredient in the determination of the primary energy of ultra-high-energy cosmic rays, 
which generate extensive air showers in the Earth's atmosphere, is a detailed model of hadronic 
multiparticle production which is reliable up to very high centre-of-mass energies, more than an 
order-of- magnitude higher than that available at the LHC. This enables the development of the 
cascade to be modelled from the primary interaction high in the atmosphere to the observations 
on the Earth. 

A recent paper [42] discusses the attributes of the Monte Carlos that have so far been used 
to generate these extensive air showers. The discussions in [42] show that none of these Monte 
Carlos are able to completely describe all the available measurements made at the high energy 
Tevatron and LHC colliders. The model we have described above has the potential to have all 
the required attributes, but it is still basically at the parton level and needs hadronization to be 
fully implemented. That is, a Monte Carlo realisation is necessary^^. However some comments 
can already be made. Above, we have discussed the multi-Pomeron structure of the QGSJETH 
Monte Carlo of Ostapchenko [35] , which has been used to generate the extensive cosmic ray air 
showers. We saw that set (A) of QGSJET predicts that (Jtot is about 30 mb larger, and, even 
more importantly, that cigjf'^^ is more than twice smaller, than the expectations of our model 
at a/s = 14 TcV. It is clear that the difference between the predictions of the two models is 
large. Moreover, for the generation of cosmic ray air showers, proton diffractive dissociation is 
very important since it produces the leading hadrons which carry a large fraction, a;^, of the 
initial energy. So confirmation, at the LHC, of the expectations for ^da/d^, and especially the 
study of low-mass dissociation, would be extremely instructive. 

7 Conclusions 

At first sight, it is surprising that the present 'soft' LHC data are tolerably well described by 
Monte Carlos based on LO DGLAP evolution. Actually, at such large collider energies we expect 
BFKL-like evolution, where there is no strong kt ordering, to be more relevant. Moreover, in 
DGLAP-based Monte Carlos, the 'soft' data are reproduced within a 'hard' approach with a 
rather high infrared cutoff k^am which increases with collider energy reaching a value /Cmin — 3 
GeV at ^ = 7 TeV. 

^^The construction of such a Monte Carlo is, at present, underway [43]. 
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In this paper we explain this puzzle. We use a perturbative QCD framework based on 
Reggeon Field Theory. That is, an approach based on the Reggc-likc structure of the BFKL 
framework in which the symmetry of the kernel allows evolution with both strong kt and 
inverse kt ordering; in other words diffusion in Inkt. We demonstrate that the rescattering and 
absorption of the intermediate partons, described by the enhanced reggeon diagrams, strongly 
suppresses the evolution in the direction of low kt, where the absorptive cross section a^^^ (xl/k^ 
is large. This acts as an effective infrared cutoff, /cmin, whose value grows with increasing energy. 
It thus explains the dominance of evolution in the 'DGLAP direction', that is with kt increasing 
from the proton to the central rapidity region. Thus we have the opportunity to describe 'hard' 
and 'soft' processes in a unified way. 

The model is based on simplified LO perturbative QCD expressions, with 'renormaliza- 
tion' parameters to account for the higher ccg order corrections, and embodies multi-Pomeron 
absorptive effects. It can be extended into the low kt domain, and it allows a satisfactory 
simultaneous description of 'soft' and 'hard' phenomena. In particular, after the parameters 
are fixed to describe the behaviour and the absolute values of the available soft cross sections 
{(Ttot,d<^ci/dt,dasY)/dM'^), the model generates quite reasonable \ow-x gluon distributions (for 
both the diffractive and inclusive processes) without fitting any additional parameters respon- 
sible for the 'input' parton distributions. Recall that, on the contrary, in conventional global 
parton analyses the 'input' distributions coming from the soft (confinement) region are com- 
pletely unknown, and their forms are parametrised and fitted to describe the experimental 
results for DIS and related 'hard' processes. 

Some predictions of the model for the LHC energies are given in Figs. 6—9 and Tables 2—5. 

Appendices 

A Evolution equations 

Here we give more details of the model that we have outhned in Section 3. We start with the 
evolution equation, (13), which generates the ladder structure of the bare Pomeron exchange 
amplitude F{y, kt, b), as shown in Fig. 10(a). The amplitude depends on the positions, ba and 
bf,, of the gluons (a, b) in the impact parameter plane, and on the available rapidity interval, 
y. It is convenient to introduce the conjugate transverse momentum variable kt to {ba — b^). 
Then the evolution equation of the amplitude in rapidity, y, may be written in the symbolic 
form 



where b = {ba + &b)/2 and K is the known BFKL-like kernel. Note that dy = d\n{l/x). 
The iteration of this equation generates the gluon ladder. This is illustrated in Fig. 10(b) for 




(35) 
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Figure 10: (a) The bare Pomeron exchange amplitude F{y, kt, b) in terms of colourless two-gluon 
exchange, where b = (b^ + bb)/2 and kt is conjugate to (6^ — b^); (b) the evolution of the bare 
Pomeron amplitude in rapidity, for the example of a proton interacting with a colourless qq— dipole 

the example of a Pomeron-exchange interaction of a proton with a quark-antiquark dipole 
positioned at ba and bb in the impact parameter plane. 
If we were to take 



then the solution is F{y) oc exp(Ay), with kf and b frozen during the evolution. Thus, this 
simple approximation gives F(Y) ~ e'^^ ~ s^. In this case the evolution equation has the 
form dF/dy = AF, where A is the probability to emit new intermediate partons per unit of 
rapidity; it is analogous to the splitting function of DGLAP evolution. 

However, the important new feature, which we explore in this paper, is the diffusion, or 
random walk, of kt as we evolve along the ladder. This is necessary in order to have a unified 
description of soft data together with, for example, high pt jet production. To allow for the 
Inkt diffusion we use a simplified form of the BFKL kernel. 



which takes into account the main qualitative features of the BFKL approach. It allows for 
diffusion in Inkt space with a diffusion coefficient controled by the parameter d which was 
chosen to reproduce the diffusion caused by the BFKL kernel with resummed next-to-leading 
ln(l/x) (NLL) corrections [24]. Here, the normalization factor, A^, is chosen so that 1 -|- A is 
the usual bare Pomeron intercept. Note that BKFL evolution is symmetric in both directions 
of rapidity, or, equivalently, is symmetric in kt and k't, when the amplitude is written in terms 
of the corresponding symmetric function [44] 



K = A 5{kt - k[)nkl 



(36) 




(37) 




(38) 
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Thus the square root factor in (37) is needed, since we consider evolution for F{k'^^) oc 
which is clearly not symmetric in kt and k'^. Recall that the exponential factor in (37) is 
symmetric. Note also that in the collinear limit k'^ <C kt (or k[ ^ kt) the kernel (37) has the 
conventional LO logarithmic behaviour K oc k[^/kf. 

At this point, it is worth emphasizing, that in the model we arc proposing, wc do not 
use the precise LO BFKL result (it is known that the NLL corrections are large), but rather 
the simplified form, (37), which takes into account the main qualitative features of the BFKL 
approach. It includes diffusion in Inkt space with a diffusion coefficient d given in an approximate 
kernel consistent with resummed next-to-leading ln(l/x) (NLL) corrections. Moreover, the 
intercept A is tuned to describe the available 'soft' data, and has a value close to that given 
by the NLL resummed BFKL equation, A ~ 0.3, see Table 1. 

Now we must include absorptive effects. That is, to allow for the contributions from multi- 
Pomeron diagrams. Since the intermediate parton may be absorbed by an interaction with the 
particles (partons) from the wave function of both the beam or target hadron, we now need to 
solve the two cquations^^ 

r)F, (ii h\ r rl'^k' 

' - ' ' exp(-A[Q,(y, fej) + Q,(y', k't)]/2) K{kt, K) F^iy, K, b). (39) 

- J ^^,2 exp(-A[Q,(y', k't) + Qk{y, Km K{kt, k't) Fi{y', k't, b). (40) 

We evolve (39) from y = yQ to y = — yo, and (40) from y' = y^ to y' = — yo; where 
y' = Yk-y and Yfc = ln(s/A;f), recall (18). 

To calculate F we need to solve these integro-differential evolution equations as a functions 
of both y and kt for a whole range of fixed values of b. This is done iteratively, giving the 
solution Fik, which depends on two indices, that is on the properties of the 'beam' and 'target' 
diffractive eigenstates. Moreover, note that the F^k now depend on two vectors in impact 
parameter space - the separation bi between the position of the intermediate parton c and the 
beam hadron, and the separation ^2 between c and the target hadron, see Fig. 11. Recall that 
the impact parameter of c is that of the 'centre', (&„ + bi,)/2, of the Pomeron. The argument 
6 in (39,40) now symbolically denotes both bi and 62- 

We have written (39) and (40) in terms of the 'simpler' form, (10), of the absorptive factors, 
in order to make the evolution equations easier to read. If form (12) is used (which is the one 
adopted in the present model), then the absorptive factor exp(— A(i7fc(?/, K) + ^i{y', fcJ))/2) in 
(39) should be replaced by the product 

xnkiy,K)/2 xn,iy',k',)/2 ■ ^ ^ 



dy J Tik't^ 

dF,{y',kt,b) f£k[ 



^^The data in the triple- Rcggcon domain indicate very small t-slopes of all the triple- Reggeon vertices [16, 45]. 
Indeed, the slopes are consistent with zero. Thus, the size of the multi-Reggeon vertices are negligible in 
comparison with the size of the incoming hadron. For this reason the absorptive corrections (that is, the 
exponential factors on the right-hand-side of (39) and (40) such that the opacities fij, Clk are taken at the same 
point in the impact parameter plane b. 
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Figure 11: The positions of the target i, beam k and the intermediate parton c in the impact 
parameter plane. 



y 



k 



yo 



Figure 12: The irreducible amplitude F^) of a high energy interaction 



B Expressions for the observables 

Here we summarize how the various amphtudes and cross sections can be calculated in terms 
of the model proposed above. 



B.l Elastic amplitude and low-mass diffraction 

To calculate the elastic amplitude we need the s-channel two-particle irreducible amplitudes for 
the scattering of the various diffractive eigenstates, i and k, for given separations 6 = 61 — 62 
between the incoming protons. These amplitudes, or effective opacities, are given by 

(6, Yk) = j d\tY j Frkiy, K K b2) F,k{Yk - y, K &i, ^2) d%id%25^''\hr -b2-b) (42) 

where Yk = ln{s/kf), see Fig. 12. Note that there is no integraP^ over y. The convolution may 
be calculated at any rapidity y, leading to the same result. Given this effective Hk eikonal', we 
can calculate the cross sections (using the analogous relation to (14)). Thus, the total cross 
section is given by 

2j2W.\'Wk\' I (1 - rf^fo , (43) 



C^tot 
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The integral over y gives the muhiphcity. 
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where Cl is shorthand for the effective opacity, (42), that is, Q = Q^^(Y,b). The a, are the 
probabihty amplitudes for the various diffractive eigenstates. That is, for the incoming beam 
proton we have \p) = X^ai|i), and similarly for the incoming target proton. Similarly for the 
elastic interaction we have 



"df ~ 47r 



i.k 



(44) 



(Tel = / d% ^^\ai\^\ak\^ [l - e 

i,k 

Allowing for the low-mass dissociation of one proton we get 



{1/2 



(45) 



Cel+SD 



d% \ai 



1 - e-^l^ 



(46) 



B.2 High- mass dissociation 

High- mass diffractive dissociation of the beam particle can be written as the elastic scattering of 
an intermediate parton c due to its absorption on the target, which is described by an amplitude 
like (14) but with O replaced by Wtkik't)- Thus, at each impact parameter point 6, and for each 
value of A;^, the cross section for single dissociation is proportional (i) to the elastic c — k cross 
section (1 — exp(— Afi^d/, 6)/2))^; (ii) to the probability to find the intermediate parton c in 
the interval dy, that is Aexp(— AQj/2 — Af2fc/2); (iii) to the amplitude Fj of the parton c-beam 
interaction; (iv) to the gap survival factor S'^{h) — exp(— Q(y, b)). The resulting cross section 
for high-mass dissociation reads 



dy 



N 



-Anfc(2/,fc^,6i,62)/2 



)Ae- 



- AHi (Ffe ,6i ,62 )/2- Af2fe (y.fc^ ,61 ,62 )/2 



X 



X -^kiy, K bi: b2)F,{Y, - y, 61, 62)^^(161 - b^Ddk.'d'hd'h , 



(47) 



where 61 (62) are the coordinates in the impact parameter plane with respect to the beam 
(target) hadron. The normalisation factor N is specified in (49). The eikonal gap survival 
probability^^ 

Slib) = eM-^fkib)) . (48) 



^^Strictly speaking, when calculating the gap survival probability in each particular case, we only have to 
account for the possibility of rescattering which produces secondaries within the gap interval. That is, in (48) 
we should not put the whole irreducible amplitude 17^^(5), but, instead, part of it; since the contribution from 
the processes with a gap in the same (or a larger) rapidity interval does not change qualitatively the structure of 
the diffractive dissociation event. In the present computations we neglect this effect. This means that actually 
the gap survival probabilities, and the true cross sections of diffractive dissociation, should be a bit larger. 
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First, we consider hypothesis (9) for the multi-Pomeron vertices g^^. Accounting for the 
different Good- Walker eigenstates, and integrating over the transverse momentum k[ of parton 
c, we obtain 



M^dasD 



^ \ak\^^JTik{y, k't, bi, b2)pik(y, K, bi, b2))Sik(\bi - 62I) 



X 



1 



X ^Vlkiy, k[, 61, b2)Fi{Yk - y, k[, 61, b2)d%id% , 
where the parton density 

_ /s.p-\ink{yM,b)+niiYk-y,k'^,b))/2 

and the elastic c — k amplitude 



(49) 
(50) 

(51) 



In an analogous way, the expressions for the i-slope of high-mass diffractive dissociation, 
for the cross section of Central Exclusive Production etc., can be obtained from corresponding 
formulae given in [21]. 

For hypothesis (11) for the multi-Pomeron vertices g^^, we find use of the AGK rules [12] 
shows that expression (49) takes the form 



M^dasD 
dM^ 



with 



\ak\'^Tik{y, K, 61, b2)\/pik{y, KMM))Sik{\bi - 62I) 



X 



X Fi{Yk - y, 61, 62, Yk - y)d%d% , 



-Pik = y(l - e-^"«(^-^''=-^))/2)/(Af7,(n - yX.^))/2). 



(52) 
(53) 



B.3 Hadron pt spectra 

From the AGK cutting rules [12], it follows that the inclusive (gluon) jet distribution is given 
just by the irreducible amplitude 



da f 71 dk 

-f-pi = \o^i?W? / ^ 2 Pikiy^ h)F^^{y, kt, 61, 62) F^^{Yk - y, kt, 61, 62) d\d%. 

y ik -'1* ^* 

(54) 

Here, following the AGK rules [12], we use the gluon densities, F_ of (26), calculated without the 
absorption caused by the diagrams which cross the rapidity of the emitted (inclusive) gluon^^; 
and account for the fact that any BFKL ladder, as in Fig. 4, with a transverse momentum 
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That is, without the absorptive effects in target opacity which arise from the beam side. 
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of the t-channel gluon k[ > qt may emit a 'soft' gluon qt with the same vertex squared factor 
|yp oc 1/qt- Therefore (54) contains an integral over k'^. 

To obtain the hadron pt spectra this distribution should be convoluted with the gluon 
fragmentation function Dg{z, qt), where z — pt/qt- We assume a simplified gluon fragmentation 
function 

Dg{z,qt) = — •(!- zf+^ with S = / Ncas{k')d]nk'yn (55) 

where the power 6 accounts for the main (double-logarithmic) effect of the scaling violation in 
the large z region. The normalization factor A'^ is fixed by energy- momentum conservation. In 
the relevant region of z this function is in satisfactory agreement with the available data on 
gluon jet fragmentation [46] . 

The results shown in Fig. 6 do not account for the secondaries coming from the hadronization 
of (the LUND or other hadronization) colour strings with small qt- Thus the deficiency of the 
prediction at low pt, seen in Fig. 6, may be removed by the particles produced via hadronization. 
To reproduce this contribution we need an explicit Monte Carlo realization of the model. 
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